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ABSTRACT 



The Schwinger-Dyson equation for the quark self-energy is solved for the case of 
■ the third and fourth quark generations. The exchanges of standard model gluons 

■^j- , and Higgs bosons arc taken into account. It is found that Higgs boson exchange 

ON ■ dominates the quark self-energy in the ultraviolet region for sufficiently large 

input quark masses causing the running quark propagator mass to increase with 
^— ^ energy-scale. No running of the quark mixing angles is found for input quark 

| masses up to and including 500 GeV. 

Unlike the renormalization group equation, the Schwinger-Dyson equation (SDE) 
^ ' analysis presented here enables us in principle to calculate the dynamical quark 

i | mass functions for very light quarks at low energies where quantum chromodynamics 

^1 (QCD) is in the non-perturbative region. This is of interest because experimentally 

accessible quantities are presently found at low energy-scales. In the presence of heavy 
quarks, such as bottom and top, it allows the calculation of running mass matrix at 
^ ■ momentum scale of the order or less than the quark masses, at which scale the RG 

equation results may not be reliable. Moreover, this nonperturbative SDE approach 
is also interesting in that it can describe the running mass function of an ultra-heavy 
quark, with a large Yukawa coupling that cannot be treated reliably by perturbation 
theory. This would be relevant for a fourth generation quark. 

In this study we assume that chiral symmetry is broken and solve the quark 
propagator SDE to find running quark mass functions and running quark mixing 
angles. In addition to the usual gluon contributions, we include Higgs boson exchange 
effects in the quark self-energy. The contributions to the quark self-energy due to 
Higgs bosons are being studied here for the first time in the context of the quark 
propagator SDE. We also study for the first time in a Higgs-boson-plus-gluon model 
the multigenerational cases of quarks which enable us to analyze quark mixing angles. 



The general form of the quark propagator SDE that we use is 

S-\q) = S \ q ) - j —g slil S{k)g s A v G^ - J g H S{k)g Y h H P H , (1) 
*Conference speaker 



where g s ^^ and g s A u are quark-gluon vertices, S(k) is the quark propagator, G^ v is 
the gluon propagator, gyA H is the quark-Higgs boson vertex factor, and Ph is the 
Higgs boson propagator. 

The ladder approximation quark propagator in the Landau gauge in Minkowski 
space that we use is iS(k) = (1/ + M){k 2 — M 2 ) -1 , where M is a 2 x 2 quark mass 
function matrix, assumed to be symmetric for this study. The gluon propagator we 
use is 

£GT{k) = i (-<T + ^pj G(k), with G(k) = ^ (2) 
and the leading-log (one-loop) coupling 

Aix 2 d 

9s(k) = ; ; », ; (3) 



ln(x — 



QCD 



where d = 12/(33 — 2rif), nf is the number of flavors. This leading- log QCD coupling 
is renormalization-scheme and gauge independent, as is its second order, two-loop, 
extension @. xq is a parameter that functions as a smooth infrared cutoff. Eq. (|3]) does 
not accurately model the gluon potential in the infrared region, but it does enable 
us to assess the relative magnitude and shape effects of the Higgs boson and gluon 
exchanges in the quark self-energy i. 

For the Yukawa sector we consider a 2 x 2 symmetric Yukawa coupling matrix, 
which can be diagonalized by a 2 x 2 unitary matrix, to give diagonal elements g±. The 
one- loop Renormalization Group Equation for the diagonalized Yukawa couplings, g±, 
without weak interaction effects, is 

W -i(^W- 2 ^ D W 9± (')) (4) 

with t = ln(xo + q 2 / Aq CD )/2, and Aqcd = 0.18 GeV. This has solution 

9± = n . + 2d 1 9 — ( 5 ) 



C±t 2d - 



16tt 2 l-2d 



where the C± must be determined by boundary conditions. We apply the boundary 
conditions m 2 ± (l GeV) = g±(l GeV)v 2 /2 for small (< 150 GeV) running quark 
masses, and m 2 ±(m±) = g±(m±)v 2 /2 for large running quark masses, where v ~ 246 
GeV. 

We solve the resulting SDE by converting it into a second order differential equa- 
tion and applying a fourth-order Runge-Kutta subroutine i. We input M bottom (0) =4.89 
GeV, M top (0) =179 GeV, and for the degenerate fourth generation quarks M bottom i (0) = 
M top '(0) =506 GeV Bi. For the third-fourth two-quark-generation case we input the 
following q 2 1 K 2 q CD = values of the mixing angles for the up-sector and down-sector: 
9 U p = 0.5 radians, and 9 down = 9 up - 9 cabi bbo with 9 cabibbo = 0.22 radians I. We are 
primarily interested in the question of the running of the mixing angles. 



Table 1. Two-quark-generation mass results, for A = 5 x 10 5 , B — 9.77 x 10 8 . 



quark 


M q (A) 


M q {B) 


Mqcd+B 
Mqcd 


bottom 


27.2A QCD 


15.3Aqcd 


1 


top 


993A QCD 


840Aq C d 


1.01 


bottom'=top' 


28WA QCD 


3150Aq C d 


7^1, varies 



The QCD plus Higgs boson case of the bottom quark mass function does not 
deviate significantly from the pure QCD running mass function. The QCD plus 
Higgs boson case of the top quark mass function shows some deviation from pure 
QCD. The Higgs boson term has an opposite sign relative to the gluon term, so it 
drives the top mass function up, in contrast to the usual QCD-only decreasing mass 
function. The Higgs boson has a substantial effect on the fourth generation quark 
self-energy. In fact, it causes the running fourth generation quark masses to increase 
in the asymptotic region. The results are summarized in the table, where A and B 
are the initial and final integration points. 

By studying the mixing angle results, we determine that the running of the mixing 
angles due to any effects in our model is less than one part in 10 9 , for any quark mass 
input value up to and including 500 GeV. This is expected for pure QCD and can 
be proven analytically for QCD i. We have been unable to prove analytically that 
the cabibbo angle is constant when the Yukawa interaction is included in the SDE 
analysis, but we have established numerically that there can be no observable effect 
even for the heavy hypothetical fourth generation. 
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